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摘要 : [H]AA AFEK Bactrocera 昆虫 是 最 具 经 济 重要 性 的 实 晶 类 害虫 ,本 研究 依据 mDNA CO 基因 碱 基 
六 列 对 离 腹 朝 毛 实 蝇 属 常见 实 晶 种 类 进行 识别 鉴定 与 系统 发 育 分 析 。【 方 法 】 以 口岸 经 常 截获 的 离 腹 寡 毛 实 蝇 属 
8 MER 21 种 实 蝇 为 对 象 ,采用 DNA 条 形 码 技术 ,通过 对 mDNA COI 基因 片段 ( 约 650 bp) 的 测序 和 比 对 ,以 
MEGA 软件 的 K2-P 双 参 数 模型 计算 种 内 及 种 间 遗 传 距离 ,以 邻接 法 (NJ) 构建 系统 发 育 树 。[ 结果 ] 聚 类 分 析 与 ] 
态 学 鉴定 结果 一 致 , 除 11 种 单一 序列 实 蝇 外 ,其 他 10 种 实 蝇 均 各 自 形成 一 个 单 系 ,节点 支持 率 为 99% 以上。 种 
(10 种 ) 遗传 距 离 为 0. 0003 ~0.0068 ,平均 为 0.0043; 种 间 (21 种 ) 遗传 距离 为 0.0154 ~ 0. 2395 ,平均 为 0.1540; 
间 遗 传 距离 为 种 内 遗传 距离 的 35.8 倍 ,而 且 种 内 、 种 间 遗 传 距离 没有 重合 区 域 。[【 结论 基于 mDNA CO 基因 
DNA 条 形 码 技术 可 以 用 于 离 腹 寡 毛 实 临 属 昆 虫 的 快速 鉴定 识别 ,该 技术 体系 的 建立 对 实 蝇 类 害虫 的 检测 监测 具 
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Abstract: [ Aim] Bactrocera is the most economically significant genus of tephritid flies. The present 
study aimed to make molecular identification and phylogeny analysis of common species of the genus 
Bactrocera. [ Methods] Twenty-one species belonging to eight subgenera of the genus Bactrocera that were 
frequently intercepted from customs were identified using DNA barcoding technology based on 
mitochondrial cytochrome c coxidase subunit I ( COI) gene. The partial fragments ( about 650 bp) of the 
mtDNA COI gene were amplified using DNA barcoding universal primers. The obtained COI gene 
fragments were sequenced and aligned. The phylogenetic tree was established by a neighbor-joining ( NJ) 
method. The intra- and inter-species genetic distances were calculated with MEGA version 5. 0 software 
using the Kimura 2-parameter model. [Results] The identification results based on the phylogenetic tree 
were consistent with those based on morphological analyses for the 21 fruit fly species. In these species, 
11 species of Bactrocera formed monophylies of themselves, the other 10 species shared one monophyly, and 
all bootstrap values of the original divergence among different haplotypes within the same species were over 
9996. The average inter-species genetic distance between the 21 species was 35. 8 times higher than the 
average intra-species genetic distance of the 10 species (0. 1540 vs. 0.0043). There was no overlap between 
intra- and inter-species genetic distances. [Conclusion] The results indicated that the DNA barcoding based 
on the partial sequence of mtDNA COI gene cauld provide rapid and accurate identification of Bactrocera 
species. The technology could be used in identifying and monitoring tephritid fruit fly species. 


Key words: Bactrocera; DNA barcoding; mtDNA COI gene; species identification ; molecular phylogeny 
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Sc Wa (fruit. fly ) HIH H ( Diptera ) Sc Wi fl 
(Tephritidae) 昆 虫 的 统称 ,是 双 翅 目 昆 虫 中 种 类 最 
多 .形态 变化 最 大 的 类 群 ,几乎 遍布 世界 各 地 ,寄主 
范围 广 , 绝 大 多 数 种 类 是 果蔬 作物 害虫 (Norbom et 
al., 2000) ,使 果蔬 产业 遭受 了 严重 的 经 济 损失 ,是 
国际 公认 的 危险 性 害虫 ,许多 国家 将 其 列 为 重要 检 



































(Walker) Wta B. (Z.) scutellaris ( Bezzi) ,.H. 
条 实则 B. (Z.) scutellata (Hendel) 、 瓜 实 量 、 桔 小 
实 蝇 和 番 石 榴 实 晶 B. (B. ) correcta (Bezzi) 等 6 种 
实 蝇 的 指纹 图 谱 ;Chua 等 (2010) 采 用 基于 COI 基因 
的 ”PCR-RFLP 
polymorphism ) 技术 成 功 地 区 分 了 杨桃 实 蝇 和 木瓜 


( restriction fragment length 








疫 对 象 。 纵 观 世 界 上 有 具 经 济 重 要 性 的 实则 主要 隶属 
5 AJE, BI FEE SZ WR JS Anastrepha Schiner , Ej i 3 -E 3E 
taj Bactrocera Macqart, € íi 3€ E X W& JE Dacus 








实 蝇 。 在 系统 发 育 分 析 方 面 ,不 少 学 者 采用 不 同 的 
分 子 标 记 技 术 对 实 晶 科 科 下 阶 元 进行 系统 发 育 关 系 
人 研究 ,如 Muraji 和 Nakahara (2001) 将 16S rRNA 分 





Fabricius | jë 3X Wg J& Ceratitis Macleay 42 3c hg Jg 
Rhagoletis Loew, HEF , 5 E I E Sce my rn 234 50 
种 实 晶 是 重要 经 济 害虫 (White and Elson-Harris, 
1992; 汪 兴 鉴 ， 1995 ) , 主要 分 布 在 亚洲 和 南 太 平 洋 
地 区 , 而 桔 小 实 晶 Bactrocera ( Bactrocera) dorsalis 
(Hendel) KERAP, Ak B. (B.) zonata 
( Saunders) iite B. (Daculus) oleae ( Gmelin) 
RJE B. (Zeugodacus) cucurbitae ( Coquillett) 
等 是 该 属 最 具 代 表 的 种 类 ( 梁 广 勤 等 , 2008), HE 
1835 F ES XE -E Sc ER Jes BL SZ 80 JP 37 OR GE , FH. ES 
1989 年 才 被 正式 确定 为 属 (Drew, 1989 ) , 且 时 至 今 
日 对 其 分 类 地 位 的 探讨 仍然 是 实 蝇 科 系统 发 育 研究 
的 重要 内 容 之 一 。 鉴 于 不 同学 者 对 该 属实 蝇 形 态 学 
的 分 类 依据 不 同 , 且 某 些 鉴 定 特征 仍 值得 商检, 进而 
在 主观 上 造成 了 该 属 分 类 的 难以 统一 (Drew，1989 ; 





















































子 标记 技术 与 12S rRNA 分 子 标记 技术 相 结合 ,分 
Pri BE A EX jE Bactrocera, Zeugodacus, 
Paratridacu 和 Gymnodacus 等 4 个 亚 属 18 种 实 晶 的 
线粒体 DNA 序列 ,从 分 子 水 平 初步 探讨 了 其 系统 发 
育 关系 ;Smith 等 (2003 ) 以 该 属 Austrodacus , Bactrocera , 
Zeugodacus , Hemigymnodacus , Daculus 和 Notodacus 等 
6 个 亚 属 25 种 实 晶 为 靶 标 ,采用 16S rRNA 和 COIT 基 
因 分 子 标记 方法 ,进行 了 物种 亲缘 关系 分 析 。Zhang 
4 (2010) 3 Fl 16S rDNA 与 COI 基因 3 ' 端 分 子 标记 
相 结合 的 方法 ,对 分 布 于 重庆 地 区 的 离 腹 寡 毛 实 蝇 
属 的 7 个 种 和 NCBI 数据 库 中 调用 的 20 个 种 的 相 
应 序列 进行 了 简单 的 系统 发 育 分 析 。 然 而 ,上 述 技 
术 方 法 或 结果 不 稳定 (如 RAPD JR) ,或 需要 酶 切 
如 RFLP 技术 ,或 缺乏 统一 的 标准 (如 16S rRNA, 
12S rRNA ,COII 等 技术 ) ,无 法 全 面 满 足 快速 ,准确 、 


























Drew and Hancock, 2000)。 此 外 ,该 属 还 拥有 众多 
在 形态 上 极其 相似 、 被 称 为 复合 种 或 近 缘 种 的 类 群 ， 
如 村 小 实 蝇 复 合 种 B. dorsalis complex 由 75 fp 3c 
组 成 (Drew and Hancock, 1994; Clarke et al., 
2005) ,长 期 以 来 被 作为 同一 物种 ,直到 后 来 才 被 陆 
续 描 述 和 定名 ;而 且 复合 种 内 各 物种 间 的 形态 特征 
非常 相似 ,种 间 特 异性 识别 特征 如 翅 脉 . 足 、 胸 部 和 
腹部 等 的 差异 相对 较 小 ,有 些 种 甚至 还 需 借助 寄主 
与 分 布 范围 等 进行 区 分 ( 陈 鹏 和 叶 辉 ,2007 ) Drew 
和 Hancock (2000) 基于 形态 学 将 离 腹 寡 毛 实 蝇 属 
Hy 28 个 亚 属 归纳 为 4 个 类 群 , 即 Bactrocera , 
Queenslandacus , Zeugodacua 和 Melanodacus > 

长 期 以 来 ,国内 外 对 实 蝇 的 研究 主要 集中 在 形 
态 学 鉴定 .生物 学 和 生态 学 特性 研究 以 及 综合 防 控 
等 方面 ( 汪 兴 鉴 ， 1995; Jamnongluk et al., 2003 ) 。 
随 着 分 子 生 物 学 技术 的 迅速 发 展 , 及 其 所 具有 的 不 
受 虫 态 或 龄 期 以 及 标本 完整 性 限制 等 特性 ,分子 标 
记 技 术 正 逐步 应 用 于 实 蝇 的 检测 鉴定 以 及 系统 发 育 
关系 人 研究。 在 检测 鉴定 方面 , 张 亮 和 张 智 英 (2007 ) 
过 RAPD 技术 构建 了 南瓜 实则 B. (Z.) tau 
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高 通 量 KUNA ES d FE fin E PEE S8 SER TR, 
此 难以 标准 化 和 规模 化 。 

基于 线粒体 DNA 细胞 色素 C 氧化 酶 亚 基 工 基 
(mtDNA COT) 的 DNA 条 形 码 ( DNA barcoding) 技 
术 即 通过 对 COI 基因 5' 端 一 段 长 度 约 为 650 bp 的 
标准 基因 片段 进行 测序 和 序列 分 析 , 进 而 在 DNA 水 
平 上 成 功 地 区 分 物种 (Hebert et al., 2003) 的 技术 ， 
其 主要 作用 包括 对 物种 进行 准确 的 识别 鉴定 ,发 现 
新 种 及 隐 存 种 ,以 及 重建 真正 的 物种 与 高 级 阶 元 的 
演化 关系 等 ,已 成 为 生物 分 类 学 研究 的 新 方向 ,对 物 
种 分 类 学 及 系统 发 育 .保护 生物 学 .生物 多 样 性 研究 
等 的 发 展 具 有 至 关 重 要 的 作用 (Valentini et al., 
2008; Floyd et ol.，2010 ) 。 目 前 ,尽管 已 有 部 分 分 
子 标记 技术 用 于 实 晶 科 昆虫 的 识别 研究 ,但 由 于 所 
用 标记 基因 的 不 同 ,其 分 析 方 法 和 结果 展示 无 法 统 
一 或 标准 化 ,因此 所 获得 的 序列 难以 实现 大 范围 共 
享 。 而 基于 COI 基因 的 DNA 条 形 码 技术 具有 统一 
性 和 标准 性 ,目前 已 成 为 生物 分 类 学 研究 中 应 用 最 
为 广泛 的 分 子 手段 ( Monaghan et al., 2005; 
Hajibabaei et al., 2006; Assefa et al., 2007 ; Shufran 
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and Puterka, 2011) ,已 用 于 鳞 翅 目 、 膜 翅 目 、 同 翅 目 
等 的 物种 鉴定 及 其 系统 发 育 关系 分 析 , 并 对 基于 形 
态 学 的 物种 分 类 进行 了 验证 或 修订 (von Dohlen et 
al., 2006; Acs et al., 2010 ; Pauls et al., 2010) ,而 基 
于 此 方法 的 实 蝇 系 统 发 育 关系 的 探讨 鲜 见 报道 。 

本 研究 围绕 离 腹 寡 毛 实 蝇 属实 蝇 种 类 多 、 形 态 
相似 、 极 易 随 果蔬 的 贸易 活动 远 距离 传播 扩散 , 却 难 
以 快速 准确 鉴别 的 问题 ,以 经 常 截 获 的 8 个 亚 属 21 
种 具有 重要 经 济 意 义 的 实 晶 及 其 COIT 基因 序列 为 
WU KH DNA 条 形 码 技术 ,进行 实 蝇 类 害虫 快速 
分 子 鉴定 识别 技术 研究 ,并 探讨 其 系统 发 育 关系 。 
由 在 为 有 效 阻 截 实 蝇 类 害虫 的 传人 以 及 进一步 传播 


























扩散 ,保护 国家 贸易 信誉 和 农业 生产 安全 提供 技术 
保障 。 





1 材料 和 方法 


1.1 供 试 实 蝇 

研究 中 所 使 用 的 21 种 离 腹 寡 毛 实 晶 属 实 蝇 的 
具体 信息 如 表 1 所 示 , 其 中 桔 小 实 量 复 合 种 (包括 
吉 小 实 蝇 杨桃 实 蝇 和 木瓜 实 蝇 ) 作为 一 个 种 进行 
分 析 ( Drew and Hancock, 1994; Clarke et al., 
2005) 。 实 晶 标 本 放 入 合 有 99% 乙醇 的 离心 管 中 ， 
4% 保存 备用 。 











表 1 供 试 实 蝇 标 本 信息 
Table 1 Fruit fly species used in the present test 



























































, aam 采集 日 期 (年 -月 测序 样本 数 
编号 种 类 地 点 xu — CEDBCOD gp i 
: . Se Collecting date Number of samples 
No. Species Locality Origin Sequence ID * 
( year-month ) sequenced 
Subgenus Bactrocera 
T BE SCR s " 
a 广东 深圳 截获 
1 B. (B.) albistrigata 2011-05 PTSY. GDSZ. 1 1 
» Shenzhen, Guangdong Intercepted 
(de Meijere) 
2 云南 元 江 FSLSY. YNYXYJ.1 — 2 
Yuanjiang, Yunnan FSLSY. YNYXYJ. 6 
Tid So 云南 保山 诱捕 SO FSLSY. YNBS. 1 - 
B. (B. ) correcta (Bezzi) Baoshan, Yunnan Trapping FSLSY. YNBS. 2 
云南 元 谋 FSLSY. YNCXYM. 1 - j 
Yuanmou, Yunnan FSLSY. YNCXYM.2 
桃 实 蝇 广东 深圳 截获 
2011-05 TSY. GDSZ.1 1 
3 B. (B.) zonata (Saunders) Shenzhen, Guangdong Intercepted S Š 
jn JR R SC 广东 深圳 截获 
4 2011-05 DASSY. GDSZ. 1 1 
B. (B.) diospyri Drew Shenzhen, Guangdong Intercepted x : 
桔 小 实 晶 复合 种 
B. (B.) dorsalis complex 
人 人 
台湾 台北 JXSY.TWTB.1- 
Taibei, Taiwan 2007:07 JXSY. TWTB. 12 12 
桔 小 实 蝇 庆 武 隆 诱捕 2008-06 JXSY. CQWL. 1 - 9 
B. (B. ) dorsalis ( Hendel) Wulong, Chongqing Trapping JXSY. CQWL. 9 
VRAK pii XSY. ZJQZ. 1 - 
浙江 衢州 xen R : 
Quzhou, Zhejiang JXSY. ZJQZ. 6 
Jes 
Vo MESS "RTI 截获 YTSY. GDSZ.1 - 
B. (B. ) carambolae 2011-05 2 
Shenzhen, Guangdong Intercepted YTSY. GDSZ.2 
Drew and Hancock 
KIEZ 
木瓜 实 蝇 ~ 东 深圳 截获 MGGSY. GDSZ. 1 - 
B. (B. ) papayae Drew 2010-07 6 
Shenzhen, Guangdong Intercepted YTSY. GDSZ.6 
and Hancock 
新 喀 里 多 尼 亚 实 蝇 URRI 截获 
6 2011-05 HJJJSY. GDSZ. 1 1 
B. (B. ) psidii (Froggatt) Shenzhen, Guangdong Intercepted SGN 
EE M ERHI He 
7 REKI Kd Bak 2011-05 — XISY. GDSZ.1 1 
B. (B. ) musae (Tryon) Shenzhen, Guangdong Intercepted 
辣椒 实 “RÝR H H- 
8 BUS ad 截获 2011-05 LJSY. GDSZ. 1 1 
B. (B. ) latifrons (Hendel) Shenzhen, Guangdong Intercepted 
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续 表 1 Table 1 continued 


采集 日 期 

















































































































编号 | ID* 测序 数量 
编号 种 类 Species 地 点 Locality 来 源 Origin Collecting date 序列 测序 数量 
No. Sequence ID * Number sequenced 
( year-month ) 
Subgenus Afrodacus 
维 斯 实 “RRI H- 
扎 维 斯 实 i T KAII BG 2011-05 ZWSSY. GDSZ.1 1 
B. (A.) jarvisi (Tryon) Shenzhen, Guangdong Intercepted 
Subgenus Austrodacus 
黄瓜 果实 晶 广东 深圳 TK 
10 n 2011-05 HGGSY. GDSZ.1 1 
B. (Au.) cucumis (French) Shenzhen, Guangdong Intercepted S 8 
Subgenus Zeugodacus 
" RORI 庆 武 隆 诱捕 jo08.07 HYSY COWL. 1 - M 
B. (Z. ) diaphora (Hendel) Wulong, Chongqing Trapping HYSY. CQWL. 10 
TE 5 Sc 广东 深圳 截获 
12 2011-05 SHDSY. GDSZ. 1 1 
B. (Z.) ishigakiensis (Shiraki) Shenzhen, Guangdong Intercepted SHDSY: GDS 
13 AUR GSY. YNYXYJ. 1 1 
Yuxi, Yunnan 
Bx 
ULL GSY. YNXSBNML. 1 1 
JR 3c Mengla, Yunnan 诱捕 DOD 
B. (Z. ) cucurbitae (Coquillett) 云南 保山 Trapping GSY. YNBS. 1 - 5 
Baoshan, Yunnan GSY. YNBS.2 
广西 南宁 GSY. GXNN. 1 - 4 
Nanning, Guangxi GSY. GXNN.4 
i IER SCR 云南 大 理 诱捕 2007-01 HQSY. YNDL. 1 - 5 
B. (Z. ) scutellaris ( Bezzi) Dali, Yunnan Trapping HQSY. YNDL.5 
15 庆 武 隆 JJGSY. CQWL. 1 - 5 
H.A& zii Wulong, Chongqing 诱捕 EG JJGSY. CQWL. 7 
B. (Z. ) scutellata ( Hendel) 云南 建 水 Trapping JJGSY. YNJS. 1 - 
Jianshui, Yunnan JJGSY. YNJS. 5 
16 云南 临 沧 NYGSY. YNLC. 1 - j 
BS 3c Lincang, Yunnan 诱捕 2007.01 NYGSY. YNLC. 7 
B. (Z.) tau (Walker) ~ 西南 宁 Trapping p NYGSY. GXNN. 1 - 2 
Nanning, Guangxi NYGSY. GXNN. 5 
Subgenus Daculus 
T 橄榄 实 蜡 “ 东 深圳 截获 TET GLGSY. GDSZ. 1 - à 
B. (D. ) oleae ( Gmelin) Shenzhen, Guangdong Intercepted j GLGSY. GDSZ. 4 
Subgenus Paratridacus 
FERM ~ 东 深 圳 截获 
18 2011-05 PZSY. GDSZ.1 1 
B. (P.) expandens (Walker) Shenzhen, Guangdong Intercepted SY UG 
Subgenus Gymnodacus 
Ti Sc wa Sc " 
l 广东 深圳 截获 
19 B. (G. ) calophylli Shenzhen. G d Int ted 2011-05 YSY. GDSZ. 1 1 
( Perkins and May) Shenzhen, Guangdong ntercepte 
Subgenus Tetradacus 
庆 武 隆 JDSY. CQWL. 1 - 
20 十 大宗 2008-07 6 
Pec uen Wulong, Chongqing 诱捕 JDSY. CQWL. 6 
> .) minax 
E T i 一 
( Enderlein) 湖北 宜 都 rapping 2009-12 JDSY. HBYD. 1 z 
Yidu, Hubei JDSY. HBYD. 6 
51 ERA SCR 云南 建 水 诱捕 ST MGDSY. YNJS. 1 - P 
B. (T. ) tsuneonis ( Miyake) Jianshui, Yunnan Trapping MGDSY. YNJS. 7 
合计 Total 126 126 


“中 国外 来 入侵 物种 数据 库 Database of Invasive Alien Species in China (DIASC) (http://www. chinaias. cn). 


1.2 DNA 提取 




















HEIC 3d A DNA 提取 参照 
法 ,并 稍 加 改进 。 分 别 取 供 试 实 晶 的 成 虫 胸部 (17 


罗 晨 等 (2002 ) Jj 


10 胸 长 ) 组 织 置 于 滴 有 20 kL 裂解 缓冲 液 (50 


mmol/L Tris-HCl, 1 mmol/L EDTA, 1% SDS, 20 
mmol/L NaCl, pH 8.0) parafilm 膜 上 ,以 PCR 管 
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底部 作为 匀 浆 器 进行 充分 研磨 , 匀 浆 液 移 人 1.5 mL 
离心 管 ;然后 以 200 pL 缓冲 液 分 4 次 冲洗 匀 浆 器 ， 
合并 混 匀 ;加 入 6 pL 和 蛋白酶 K(20 mg/mL) ,充分 混 
匀 后 于 65 仿 水浴 1 h( 中 途 混 匀 3 次 ) ;加 入 220 pL 
氯仿 / 异 戊 醇 (v: v = 24: 1) ,轻柔 混 匀 后 , 冰 浴 30 
min ;然后 4% 12 000 r/min 离心 20 min, XC EY 
液 ,加 入 440 pL 预 冷 无 水 乙醇 ,轻柔 混 匀 后 于 
-40*C 放置 30 min;4%C 12 000 r/min 离心 15 min, 
小 心 弃 去 上 清 液 。 加 入 440 pL 预 冷 的 75% 乙醇 洗 
涤 ,4% 12 000 r/min 离心 15 min ,小 心 弃 去 上 清 液 ; 
然后 将 离心 管 倒 扣 于 洁净 滤纸 上 ,自然 干燥 。 每 管 
加 入 30 uL 超 纯 水 ,充分 溶解 后 于 -40% 保存 备用 。 
1.3 PCR 扩 增 .电泳 检测 及 序列 测定 

以 实 蝇 胸部 组 织 DNA. 为 模板 ,以 DNA 条 形 码 
用 型 引物 LCO1490/HCO2198 (上 游 引 物 
LCO1490 的 碱 基 序列 为 5'-GGTCAACAAATC- 
ATAAAGATATTGG-3 ,下游 引物 HCO2198 的 碱 基 
序列 为 5'-CTAAACTTCAGGGTGACCAAAAAATCA- 
3') (Folmer et al., 1994) (由 上 海 生 工 生物 工程 技术 
服务 有 限 公 司 协助 合成 ) 进行 PCR. 扩 增 ,产物 大 小 
约 为 700 bp。 反 应 体系 为 25 pL, 其 中 超 纯 水 9. 5 
uL,2 x EasyTaq PCR SuperMix 12. 5 &L,10 pmol/L 
上 游 引 物 和 下 游 引 物 各 0.5 uL, DNA 模板 2 pL, 
扩 增 程序 为 :94% 预 变性 3 min;94%C45 s,48%C 45 
s,72% 45 s,35 个 循环 ;最 后 72C 延伸 10 min, 
PCR 扩 增 反应 在 ABI-9700 基因 扩 增 仪 上 进行 , 取 4 
pL PCR 扩 增 产物 在 1. 5% (重量 ) 琼脂 糖 
(Amresco) 凝 胶 上 以 80 V 电泳 分 离 45 min ,然后 以 
Gel-Doc Universal Hood II 型 凝 胶 成 像 系统 (Bio-Rad 
Laboratories) 分 析 电 泳 结果 。 然 后 ,对 电泳 检测 合格 
的 实 蝇 DNA 的 PCR 扩 增 产物 进行 纯化 ,由 北京 三 
博 远 志 生 物 技 术 有 限 责 任 公司 协助 完成 测序 。 根 据 
标本 数量 以 及 PCR 产物 的 质量 ,每 种 实 蝇 检 测 1 ~ 
35 头 ( 表 1) ,共计 检测 126 头 , 并 将 所 有 的 126 条 序 
列 录 入 中 国 主要 外 来 人 侵 昆虫 DNA. 条 形 码 识别 系 
统 ( Database of Invasive Alien Species in China, 
DIASC) (http://www. chinaias. en) (Æ 1) ; 
1.4 COI 基因 序列 处 理 与 分 析 

以 序列 阅读 软件 Chromas 1.0 读 取 序列 并 进行 

ATHER, 以 确保 序列 的 高 质量 。 然 后 ,用 Bioedit 
7.0 软件 对 序列 进行 整合 处 理 , 以 ClustalW Multiple 
Alignment 软件 进行 对 齐 处 理 ; 并 利用 NCBI, BOLD 
等 数据 库 ,对 所 获得 的 高 质量 COT 基因 序列 与 数据 
库 中 已 公开 的 实 蝇 类 碱 基 序列 进行 同 源 性 比 对 分 
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析 。 然 后 , 以 地 中 海 实 蝇 (GenBank 登录 号 : 
DQ116368 ) 和 非洲 芒果 实 蝇 ( GenBank 登录 号 : 
AY788421 ) 为 外 群 , 利 用 MEGA 5. 0 软件 (Tamura et 
al., 2011) 以 Kimura 2-parameter ( K2-P) 双 参 数 模 型 
(Kimura, 1980) 计算 种 内 及 种 间 遗 传 距 离 ,并 采用 
邻接 法 (neighbor-joining, NJ) 构 建 系统 发 育 树 ,系统 
发 育 树 分 支 的 置信 上 度 采 用 自 展 法 (bootstrap 
analysis, BP) 重复 检测 1 000 次 。 然 后 ,结合 NCBI 
中 已 公开 的 离 腹 窒 毛 实 晶 属 35 种 重要 实 晶 类 害虫 
的 高 质量 COIT 基因 序列 ( 表 2) ,与 本 研究 的 21 种 实 
晶 ( 其 中 有 12 种 实 蝇 为 本 研究 独 享 种 类 ) 的 COI 基 
因 序 列 , 各 选取 一 条 序列 , 即 共 计 47 种 实 蝇 的 COI 
基因 序列 ,以 来 自 数据 库 中 共享 种 类 的 序列 作为 标 
准 序列 ,采用 邻接 法 构建 系统 发 育 树 ,BP 重复 检测 
1 000 次 。 

此 外 ,以 本 研究 的 3 种 桔 小 实 蝇 复合 种 为 靶 标 ， 
整合 数据 库 中 已 公开 的 其 他 4 种 桔 小 实 晶 复 合 种 
(由 于 目前 数据 库 中 尚未 有 沙 梨 果实 蝇 COIT 基因 序 
列 信息 的 收录 , 故 未 包括 其 中 ) 的 高 质量 COT 基因 
序列 ( 亦 即 片段 大 小 与 本 研究 的 靶 标 片段 基本 一 
致 )( 表 3) ,经 对 齐 处 理 后 ,进行 桔 小 实 量 复 合 种 的 
种 内 及 种 间 遗 传 距离 分 析 。 用 于 分 析 的 序列 共计 
72 条 ,其 中 村 小 实 蝇 38 条 、 木 瓜 实 蝇 10 条 ,杨桃 实 
晶 10 2&  3ERETEARSCNG 7 条 、 芒 果实 蝇 5 条 、 胡 桃 实 
蝇 1 条 、 康 迪 果 实 蝇 1 条。 然后 ,从 各 复合 种 中 随机 
选取 1 条 序列 ,以 地 中 海 实 蝇 (DQ116368 ) 为 外 群 ， 
以 邻接 法 构建 系统 发 育 树 ,BP 重复 检测 1 000 次 。 


2 结果 与 分 析 



















































































2.1 PCR 扩 增 、 序 列 测定 及 其 变异 分 析 

对 经 电泳 检测 验证 合格 的 PCR 产物 进行 纯化 
和 碱 基 序列 测定 后 ,以 Bioedit 7. 0 软件 进行 整合 处 
理 ,ClustalW Multiple Alignment 软件 进行 对 齐 处 理 ， 
修剪 成 长 度 为 648 pp 的 片段 。 然 后 ,对 所 获得 的 高 
质量 COI 基因 序列 与 数据 库 ( NCBI, BOLD ) 中 已 公 
开 的 实 晶 类 昆虫 碱 基 序 列 进行 同 源 性 比 对 分 析 , 结 
果 显 示 , 本 研究 中 21 种 实 蝇 的 COT 基因 的 碱 基 序列 
与 数据 库 中 已 知 的 训 标 种 相对 应 的 COT 基因 片段 
的 一 致 性 为 98. 90% ~ 100% , 且 不 存在 碱 基 插 和 信和 
缺失 现象 ,表明 所 获得 的 COL 基因 序列 准确 可 靠 。 
采用 MEGA 5.0 软件 ,对 21 PPk 648 bp 的 COI 
基因 的 碱 基 序 列 进行 分 析 , 结 果 表 明 , 保 守 位 点 C 


( conserved sites )406 个 ,变异 位 点 V(variablesites ) 
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表 2 用 于 DNA 条 形 码 技术 研究 的 NCBI 中 已 公开 的 35 种 离 腹 嘉 毛 实 蝇 属 实 蝇 
Table 2  Thirty-five fruit fly species of genus Bactrocera from NCBI used in DNA barcoding 






























































编导 种 类 GenBank 登录 号 “| 编号 种 类 GenBank 登录 号 
No. Species GenBank accession no. No. Species GenBank accession no. 
IR SWR Bactrocera 小 昆士兰 实 蝇 
1 DQ116237 20 < DQ116304 
(Austrodacus) cucumis ( French) Q B. (B. ) neohumeralis ( Hardy ) Q 
扎 维 斯 实 晶 
2 DQ11628 21 B. (B. igrotibialis ( Perkins DQ006870 
B. (Afrodacus) jarvisi ( Tryon) g ? UB) migrotibialis (Perkins) Q 
3 B. (Af. ) minuta (Drew) DQ116301 22 B. (B.) quadrisetosa ( Bezzi) DQ116327 
TEX BSCRSCHR 
4 DQ116234 23 B. (B. ilineola Dre DQ116338 
B. (Bactrocera) aquilonis ( May) Q (o) ihre 22x Q 
z RAR SE B. B. 3 S 
5 FEARS h (B. ) ede DQ116236 24 巴布亚 新 儿 内 实 蝇 DQ116342 
(Hardy & Adachi) B. (B. ) trivialis ( Drew) 
HRK 昆士兰 实 蝇 
6 DQ11623 25 DQ116336 
B. (B. ) cacuminata (Hering) Q 2 B. (B.) tryoni ( Froggatt) Q 
7 B. (B. ) cognata (Hardy & Adachi) DQ116251 26 EGG DQ116348 
B. (B. ) umbrosa ( Fabricius) 
dn M RES B. (B. ) verbascifoliae 
DQ116264 27 DQ11634 
3 B. (B. ) correcta ( Bezzi) g Drew & Hancock g ? 
LEMI 
9 An Dii ca DQ116261 28 ESRB. (B.) zonata (Saunders) DQ116362 
B. (B. ) curvipennis ( Froggatt) 
面包 果实 网 
10 DQ116282 2 B. (Daculu: 'anus (M GQ505008 
B. (B. ) distincta (Malloch) Q 9 ( Daculus) montyanus ( Munro) Q505 
11 tihi B. (B. ) dorsalis ( Hendel) DQ116275 30 B. (D. ) munroi White FJ009204 
ORI 
12 DQ116284 31 HN SNR B. (D.) ol ;meli DQ011886 
B. (B. ) endiandrae (Perkins & May) g MORS 人 Q 
1 JN Sc 
13 DQ116285 32 B. (Gymnodacu: l M 30504997 
B. (B. ) facialis ( Coquillett) Q ( Gymnodacus) amplexa ( Munro) GQ50499 
Tft Scu JNE B. (Zeugodacus) 
14 DQ116286 33 DQ116245 
B. (B. ) frauenfeldi (Schiner) g cucurbitae ( Coquillett ) Q 
15 ARSCH B. (B.) invadens Drew FJ009202 34 WIRE B. (Z. ) tau (Walker) DQ116245 
黄 侧 条 实 蝇 B. 
16 RÆK B. (B.) kirki (F DQ116294 35 die DQ116355 
NAROESEHR 人 Q 4 ( Notodacus) xanthodes ( Broun) g 
T SEHR. Ceratitis 
17 ”辣椒 实 晶 B. (B.) latifron (Hendel) FJ009203 36 Harp unes DQ116368 
( Ceratitis) capitata ( Wiedemann) 
B. (B. kili sis MTER C. 
T (B. ) makilingensis DQ116302 37 非洲 芒果 实 凯 AY788421 
Drew & Hancock ( Ceratalaspis) cosyra ( Walker) 
1 DQ1162 
? B. (B. ) melanotus ( Coquillett ) g 2 
* 外 群 Outgroup. 
表 3 7 种 桔 小 实 蝇 复 合 种 的 COI 基因 片段 相关 信息 
Table3 The COI gene information of seven species of Bactrocera dorsalis complex 
编号 种 类 GenBank 登录 号 来 源 
No. Species GenBank accession no. Origin 
1 胡桃 实则 Bactrocera (Bactrocera) caryeae ( Kapoor) DQ116266 NCBI 
2 TRW B. (B.) occipitalis (Bezzi) DQ116311 NCBI 
3 Wt B. (B. ) carambolae Drew and Hancock YTSY. GDSZ. 1 DIASC 
4 村 小 实则 B. (B. ) dorsalis ( Hendel) JXSY. TWTB. 1 DIASC 
5 RR B. (B. ) kandiensis Drew and Hancock DQ116295 NCBI 
6 KJK HE B. (B. ) papayae Drew and Hancock MGGSY. GDSZ. 1 DIASC 
7 JEE RK B. (B.) philippinensis Drew and Hancock DQ116317 NCBI 
8 Ji pif SES Ceratitis (Ceratitis) capitata ( Wiedemann) * DQ116368 NCBI 


“外 群 Outgroup. DIASC; 中 国外 来 人 侵 物 种 数据 库 Database of Invasive Alien Species in China (http://www. chinaias. cn). 
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240 个 ,简约 信息 位 点 Pi( parsimony information 
sites) 218 ^, FUE TT BR Ze A E I SNPs (single 
nucleotide polymorphisms ) 22 个 。 

2.2 遗传 距离 

2.2.1 离 腹 穿 毛 实 蝇 属 21 种 实 晶 的 遗传 距离 : 采 
用 软件 MEGA 5. 0 中 的 K2-P 双 参 数 模型 计算 离 腹 
宅 毛 实 蝇 属 实 蝇 的 种 内 (10 种 ) 与 种 间 (21 种 ) 遗 传 
距离 ,结果 显示 ,种 内 遗传 距离 为 0.0003 ~ 0. 0068, 
平均 为 0.0043; 种 间 遗 传 距离 为 0. 0154 ~ 0. 2395, 
平均 为 0. 1540。 其 中 ,平展 实 晶 和 艳 实 晶 间 的 遗传 
距离 最 大 ,为 0.2395 ;石垣 岛 实 蝇 和 具 条 实 蝇 之 间 
的 遗传 距离 最 小 ,为 0.0154( 表 4)。 而 且 , 种 内 遗 
传 距离 与 种 间 遗 传 距 离 没 有 重 共 区域 (图 1)。 
2.2.2 村 小 实 晶 复 合 种 的 遗传 距离 :以 本 研究 中 的 
3 种 桔 小 实 蝇 复 合 种 和 数据 库 中 已 公开 的 其 他 4 种 
桔 小 实 蝇 复 合 种 的 COI 基 因 序 列 ( 处 理 后 的 片段 
大 小 为 612 bp) WIER, XH MEGA 5.0 中 的 K2-P 
双 参 数 模 型 计算 桔 小 实 量 复 合 种 的 种 内 (5 种 ) 和 种 
间 (7 Rp) 遗传 距离 。 结 果 显 示 , 种 内 遗传 距离 为 
0. 0006 ~ 0. 0080 ,平均 为 0.0059 ,其 中 菲律宾 果实 
晶 的 种 内 遗传 距离 最 小 ,杨桃 实 蝇 的 种 内 遗传 距离 
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= 种 内 Intra-species 
= 种 间 Inter-species 


分 布 频率 Frequency (96) 
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遗传 距离 Genetic distance (x10?) 

图 1 常见 离 腹 寡 毛 实 蝇 属 实 蝇 的 种 内 (10 种 ) 与 
种 间 (21 种 ) 遗 传 距离 的 分 布 频率 


Frequency distribution of the K2-P intra-species 


8 10 12 14 16 18 20 22 














Fig. 1 
(10 species) and inter-species (21 species) genetic 


distances of fruit fly species of the genus Bactrocera 


最 大 ( 表 5)。 村 小 实 蝇 复 合 种 的 种 间 遗 传 距 离 为 
0. 0040 ~ 0. 0816 ,平均 为 0.0428; 其 中 ,杨桃 实 蝇 和 
胡桃 实 蝇 之 间 的 遗传 距离 最 大 , 非 律 宾 果 实 蝇 与 村 
小 实 蝇 及 木瓜 实 蝇 之 间 的 遗传 距离 最 小 ,而 村 小 实 
蝇 与 木瓜 实 蝇 之 间 的 遗传 距离 为 0.0067( 表 5)。 




















表 5 桔 小 实 蝇 复 合 种 的 种 内 与 种 间 遗 传 距离 


Table 5 Intra- and inter-species genetic distances of seven species of Bactrocera dorsalis complex 


种 间 Inter-species 





桔 小 实 蝇 杨桃 实 蝇 


种 内 





胡桃 实 蝇 


芒果 实 蝇 。 康 迪 果 实 蝇 。 木 瓜 实 蝇 JERR 











Intra-species 


Bactrocera dorsalis B. carambolae B. caryeae B. occipitalis B. kandiensis B. papayae B. philippinensis 

















桔 小 实 晶 Bactrocera dorsalis 0. 0068 
杨桃 实 晶 B. carambolae 0. 0112 0. 0080 
胡桃 实 晶 B. caryeae 0.0787 0.0816 - 
RKW B. occipitalis 0. 0187 0. 0256 0. 0802 0. 0078 
康 迪 果实 蝇 B. kandiensis 0. 0741 0.0715 0. 0165 0. 0737 - 
KINH B. papayae 0. 0067 0. 0116 0. 0755 0. 0184 0. 0709 0. 0063 
JEE ER B. philippinensis 0. 0040 0. 0101 0. 0776 0. 0172 0. 0730 0. 0040 0. 0006 
- : 单一 样本 ,不 适宜 进行 种 内 遗传 距离 计算 One sample only, intra-species genetic distance calculation is not available. 


2.3 系统 发 育 树 构建 
运用 MECA 5. 0 分 析 软 件 , 以 地 中 海 实 蝇 和 非 
洲 芒果 实 蝇 为 外 群 ,以 21 种 实 蝇 共计 126 条 序列 




















有 较 近 的 亲缘 关系 。 此 外 ,Zeugodacus 亚 属 为 并 系 
类 群 ,石垣 岛 实 蝇 、 具 条 实 蝇 、 黑 漆 果 实 蝇 和 黑 颜 
果实 晶 聚 为 一 大 支 , 而 南瓜 实 晶 、 瓜 实 晶 和 黄瓜 果 
























































为 靶 标 ,采用 邻接 法 构建 系统 发 育 树 。 结 果 显 示 ， 
聚 类 分 析 与 形态 学 鉴定 结果 一 致 , 除 11 种 单一 序 
列 的 实 晶 以 外 ,其 他 10 种 实 晶 均 各 自 形 成 一 个 单 
系 (图 2), 节 点 支持 率 在 99% 以 上 ,说 明 DNA 条 
形 码 可 对 离 腹 窒 毛 实 蝇 属 的 实 蝇 进 行 准确 的 识别 
和 鉴定 。 同 时 ,系统 发 育 分 析 还 显示 ,Tetradacus 
亚 属 的 柑 村 大 实 蝇 和 蜜柑 大 实 蝇 与 Paratridacus W. 
属 的 平展 实 蝇 聚 为 一 支 ( 图 2) ,表明 ,两 个 亚 属 具 
































KER N— AKE, 其 中 南瓜 实 蝇 与 瓜 实 蝇 为 姐妹 
群 (图 2)。 

从 本 研究 测 得 的 21 种 实 蝇 的 COI 基因 序列 ,以 
及 NCBI 中 已 公开 的 离 腹 寒 毛 实 蝇 属 的 35 种 重要 
实 蝇 类 害虫 的 高 质量 COT 基因 序列 ( 表 3) ,各 选 一 
条 , 即 共 10 个 亚 属 47 种 实 蝇 的 COI 基因 序列 ,以 邻 
接 法 构建 系统 发 育 树 ,结果 显示 每 种 实 蝇 均 具有 一 
独立 分 支 ( 图 3)。 
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Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 






JXSY.ZJQZ.2 
JXSY.TWTB.4 
JXSY.ZJQZ.4 


JXSY.ZJQZ.6 


JXSY.ZJQZ.1 


Bactrocera oleae |GLSY.GDGZ.1 
Bactrocera oleae | GLSY.GDGZ.4 
Bactrocera oleae | GLSY.GDGZ.3 
Bactrocera oleae |GLSY.GDGZ.2 
Bactrocera diospyri |DASSY.GDSZ.1 





Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera dorsalis 
Bactrocera papayae |MGSY.GDGZ.2 

Bactrocera dorsalis |JXSY.CQWL.2 

Bactrocera dorsalis |JXSY.COWL.6 

Bactrocera dorsalis |JXSY.COWL.4 

Bactrocera dorsalis |JXSY.COWL.7 

Bactrocera dorsalis JXSY.CQWL.3 

Bactrocera dorsalis |JXSY.TWTB.1 

Bactrocera dorsalis |JXSY.CQWL.1 

Bactrocera dorsalis |JXSY.TWTB.12 

Bactrocera dorsalis |JXSY. TWTB.11 

Bactrocera papayae |MGSY.GDGZ.6 

Bactrocera papayae |MGSY.GDGZ.3 

Bactrocera papayae |MGSY.GDGZ.4 

Bactrocera papayae |MGSY.GDGZ.1 

Bactrocera dorsalis |J XSY. TWTB.9 

Bactrocera dorsalis |JXSY. TWTB.7 

Bactrocera carambolae |YTSY.GDSZ.2 

Bactrocera carambolae | YTSY.GDSZ.1 

Bactrocera papayae |MGSY.GDGZ.5 

Bactrocera musae |XJS Y. GDSZ.1 

Bactrocera calophylli | YSY.GDSZ.1 
Bactrocera psidii |HAJJJSY.GDSZ.1 

Bactrocera albistrigata |PTSY.GDSZ.1 

Bactrocera zonata |TSY.GDSZ.1 

654 Bactrocera correcta |FSLSY.YNBS.1 

* Il Bactrocera correcta |FSLSY.YNYXYJ.5 

9L Bactrocera correcta |FSLSY. YNYXYJ.2 

Bactrocera correcta |FSLSY.YNBS.2 

Bactrocera correcta |FSLSY. YNYXYJ.3 

9| | Bactrocera correcta |FSLSY. YNYXYJ.1 

Bactrocera correcta |FSLSY. YNCXYM.2 

s| Bactrocera correcta |FSLSY. YNCXYM.1 
Bactrocera correcta | FSLSY. YNYXYJ.6 

Bactrocera correcta |FSLSY. YNYXYJ.4 

Bactrocera jarvisi |ZWSSY.GDSZ.1 

Bactrocera latifrons |LJSY.GDSZ.1 

Bactrocera minax |JDSY. HBYD.4 
Bactrocera minax |JDSY.CQWL.2 

5} Bactrocera minax |JDSY.COWL.6 
»|j Bactrocera minax |JDSY.CQWL.3 

v! Bactrocera minax |JDSY.CQWL.1 
Bactrocera minax |JDSY.HBYD.6 
Bactrocera minax |JDSY. HBYD.2 
Bactrocera minax |JDSY.HBYD.1 

3 Bactrocera minax |JDSY.HBYD.3 

3p Bactrocera minax |JDSY.CQWL.5 
Bactrocera minax |JDSY.HBYD.5 
Bactrocera minax |JDSY.CQWL.4 
€», Bactrocera tsuneonis |MGDSY.YN.5 
Bactrocera tsuneonis |MGDSY.YN.2 


JXSY.ZJQZ.5 
JXSY.ZJQZ.3 
JXSY.TWTB.10 
JXSY.TWTB.5 
JXSY.TWTB.3 
JXSY.TWTB.2 
JXSY.CQWL.9 
JXSY.TWTB.6 
JXSY.CQWL.8 
JXSY.TWTB.8 
|JXSY.CQWL.5 
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Bactrocera tau |NYGSY.GXNN.1 
Bactrocera tau |NYGSY.YNLCYX.2 
Bactrocera tau |NYGSY.GXNN.5 

1s| Bactrocera tau |NYGSY.GXNN.2 
Bactrocera tau |NYGSY. YNLCYX.3 

33] Bactrocera tau |NYGSY.YNLCGM.1 

zf Bactrocera tau |NYGSY. YNLCGM.3 
Bactrocera tau |NYGSY.GXNN.4 
Bactrocera tau |NYGSY. YNLCGM.4 
Bactrocera tau |NYGSY. YNLCYX.1 
Bactrocera tau |NYGSY.YNLCGM.2 
Bactrocera tau |NYGSY.GXNN.3 
Bactrocera cucurbitae |GSY. YNXSBNML.1 
Bactrocera cucurbitae |GSY.GXNN.2 
Bactrocera cucurbitae |GSY.YNBS.1 
Bactrocera cucurbitae |GSY. YNYXYJ.1 
Bactrocera cucurbitae |GSY.GXNN.3 
Bactrocera cucurbitae |GSY.GXNN.1 
Bactrocera cucurbitae |GSY.GXNN.4 
Bactrocera cucurbitae |GSY. YNBS.2 
Bactrocera cucumis |HGSY.GDSZ.1 
Bactrocera ishigakiensis |SYDSY.GDSZ.1 
9f Bactrocera scutellata |JTSY. YNJS.3 
Bactrocera scutellata |JTSY.YNJS.4 
Bactrocera scutellata JTSY.CQWL.5 
Bactrocera scutellata |JTSY.CQWL.7 
Bactrocera scutellata |JTSY.CQWL.2 
Bactrocera scutellata |JTSY.COWL.6 
Bactrocera scutellata |JTSY.COWL.3 
Bactrocera scutellata |JTSY.YNJS.1 

€] Bactrocera scutellata |JTSY.YNJS.2 
Bactrocera scutellata JTSY.COWL.4 
Bactrocera scutellata |JTSY.YNJS.5 
Bactrocera scutellata |JTSY.CQWL.1 
Bactrocera scutellaris |HOSY.YNDL.2 
Bactrocera scutellaris |HOSY.YNDL.4| 
71 Bactrocera scutellaris |HQSY.YNDL.5 

24 Bactrocera scutellaris |HOSY. YNDL.1 
Bactrocera scutellaris |HOSY.YNDL.3 
Bactrocera diaphora |HYSY.COWL.9 
Bactrocera diaphora |HYSY.COWL.3 
Bactrocera diaphora |HYSY.CQWL.2 
Bactrocera diaphora |HYSY.CQWL.7 
Bactrocera diaphora |HYSY.CQWL.5 
el Bactrocera diaphora |HYSY.CQWL.4 
Bactrocera diaphora |HYSY.CQWL.10 
4| Bactrocera diaphora |HYSY.CQWL.8 
Bactrocera diaphora |HYSY.CQWL.6 
Bactrocera diaphora |HYSY.CQWL.1 
Ceratitis capitata |DZHSY.GDGZ.1 

Ceratitis cosyra |FZMGSY.GDSZ.1 








Bactrocera 
99| Bactrocera 
Bactrocera 
Bactrocera 
Bactrocera 


tsuneonis |MGDSY.YN.7 
tsuneonis |MGDSY.YN.4 
tsuneonis |MGDSY.YN.1 
tsuneonis |MGDSY.YN.6 
tsuneonis |MGDSY.YN.3 






Bactrocera expandens |PZSY.GDSZ.1 








图 


Fig.2  Neighbor-joining dendrogram based on the analysis of COI gene sequences of 21 fruit fly species of the genus Bactrocera 





分 支 处 上 方 数值 表示 重复 1 000 次 后 的 自 展 值 ( > 50% ) ;标尺 示 遗 传 距离 ; 


2 邻接 法 构建 的 21 种 离 腹 寡 毛 实 晶 


属实 晶 的 系统 发 育 树 








4 [i]; Numbers on branches are bootstrap values ( > 50% ) 


图 3 和 图 








under 1 000 replicates; the scale bar represents the genetic distance. The same for Figs. 3 and 4. 


VAR pig Sc dg N PME. RRRA RDK Mg 
合 种 的 7 种 最 具有 经 济 重 要 性 的 实 蝇 的 COL 基因 序 
列 进行 系统 发 育 分 析 , 结 果 显 示 ,7 种 桔 小 实 晶 复合 
种 罕 为 两 支 ,其 中 桔 小 实 蝇 杨桃 实 蝇 \ 木 瓜 实 蝇 、 芒 
RNE .菲律宾 果实 晶 聚 为 一 支 , 胡 桃 果实 蝇 和 康 迪 
RSCHRRON 3 — 5 CELA) 。 
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备 两 个 先决 条 件 :第 一 
条 形 码 , 且 不 同 物种 间 无 共享 


讨论 


DNA 条 形 码 技术 对 物种 的 准确 识别 鉴定 需 具 
,每 个 物种 具有 独特 的 DNA 
型 序列 ;第 二 ,种 间 差 
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B.(Bactrocera) dorsalis JXSY. COWL.1 
B. (B.) dorsalis |DQO116275 
B. (B.) invadens |FJ009202 
B. (B.) verbascifoliae |DO116349 
B. (B.) cacuminata |DQ116239 
B. (B.) makilingensis |DO116302 
B. (B.) musae |XJSY.GDSZ.1 
B. (B.) zonata |TSY.GDSZ. 1 
B. (B.) zonata |DO116362 
B. (B.) correcta |FSLSY.YNYXYJ.1 
wL B.(B.) correcta |DO116264 
wor B. (A.) jarvisi |Z:WSSY.GDSZ.1 
站 L B. (Afrodacus) jarvisi| DO116289 
B. (B.) nigrotibialis |DQ006870 
Eh LF— — —— B. (B.) melanotus |DQ116299 


34 B. (B.) kirki |DO116294 
I] easet .[— B. (B.) psidii |HJJJSY.GDSZ.1 
z| —— B. (B.) trilineola |DQ1 16338 


r— B. (B.) albistrigata |PTSY.GDSZ.I1 
:L— B. (B.) frauenfeldi |DO116286 
B. (B.) arecae |DO116236 

B. (B.) distincta |DO116282 















R 








s 








E B. (B.) endriandrae |DQ116284 
4 —————— B. (Gymnodacus) calophylli |YSY.GDSZ.1 
$ B. (B.) umbrosa |DQ116348 
z B. (B.) facialis |DO116285 
HEBES B. (Af.) minuta |DQ116301 


E B. (B.) quadrisetosa |DO116327 
B. (B.) trivialis |DO116342 
B. (B.) cognata |DQ116251 
r B. (B.) curvipennis DO116261 
Em B. (B.) aquilonis |DQ116234 
s| B. (B.) neohumeralis |DO116304 
s! B.(B.) tryoni |DO116336 
r B. (B.) latifrons |FJ009203 
r 1w! B. (B.) latifron |LJSY.GDSZ.1 
B. (Gymnodacus) amplexa |GO504997 
B. (Notodacus) xanthodes |DQ116355 
B. (B.) diospyri ||DASSY. GDSZ.1 
B. (T.) minax |JJDSY.HBYD.1 
B. (T.) tsuneonis |MGDSY.YN.1 
B. (P.) expandens PZSY.GDSZ.1 
B. (Daculus) montyanus |GQ505008 
B. (D.) munroi |FJ009204 
B. (D.) oleae |DQ011886 
B. (D.) oleae | GLSY. GDGZ.1 
B. (Zeugodacus) diaphora |HYSY. COWL.10 
B. (Z.) scutellata JISY.YNJS. 1 
B. (Z.) ishigakiensis |SYDSY.GDSZ.1 
B. (Z.) scutellaris |ÁHOSY.YNDL.1 
B. (Austrodacus) cucumis |DO116237 
B. (Austrodacu) cucumis | HGSY.GDSZ.1 
w] B. (Zeugodacus) cucurbitae |DO116245 
B. (Z) cucurbitae |GSY.GXNN.1 
B. (Z.) tau |GO154160 
B. (Z) tau |NYGSY.YNLCYX.1 
[— — —— Ceratitis (Ceratitis) capitata |DO116368 
100 C. (Ceratalaspis ) cosyra |AY788421 
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图 3 ”邻接 法 构建 的 47 PiE EKR KAR RRRA W 
Fig.3  Neighbor-joining dendrogram based on the analysis of COI gene 





sequences of 47 fruit fly species of the genus Bactrocera 


















































sr Bactrocera dorsalis |JXSY.ZJQZ.1 异 应 远大 于 种 内 差异 ( Hebert et al., 2004a, 

Buc qu. PL RE SCHRUR IU] 21 种 重要 实 蝇 类 害虫 进行 识别 

ml. p bouis DO 1295 鉴定 ,并 对 其 进行 系统 发 育 分 析 , 以 K2-P 双 参 数 模 

Ceratitis capitata|DQ116368 型 计算 种 内 (10 种 ) 及 种 间 (21 种 ) 遗传 距离 。 结 果 

XU 显示 ,种 内 遗传 距离 为 0. 0003. ~ 0. 0068 ,平均 为 
图 4 邻接 法 构建 的 7 种 桔 小 实 晶 复 合 种 的 系统 发 育 树 0. 0043 ;种 间 遗 传 距离 为 0.0154 ~ 0. 2395 ,平均 为 
Fig.4 Neighbor-joining dendrogram based on 0. 1540 ; 种 间 遗 传 距离 为 种 内 遗传 距离 的 35.8 倍 ， 

the analysis of COI gene sequences of seven 完全 符合 种 间 差 异 应 大 于 种 内 差异 10 倍 的 物种 鉴 


species of Bactrocera dorsalis complex 


定 规则 (Hebert et al., 2004a, 2004b) , Aliabadian 
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等 (2009 ) 提 出 ,物种 的 准确 鉴定 与 靶 标 物种 的 种 间 
遗传 距离 与 种 内 遗传 距离 之 间 的 分 离 程度 有 关 , 当 
种 间 遗 传 距 离 和 种 内 遗传 距离 不 存在 重 辣 区 域 时 ， 
便 能 准确 地 鉴定 物种 。 尽 管 石 垣 岛 实 蝇 和 具 条 实 晶 
间 的 遗传 距离 为 0.0154 ,小 于 条 形 码 有 效 性 的 检验 
标准 , 即 种 间 遗 传 差异 应 大 于 0.02 (Hebert et al., 
2003) ,但 本 人 研究 的 种 内 与 种 间 遗 传 距离 分 布 频率 
的 分 析 结 果 显 示 ,种 内 遗传 距离 与 种 间 遗 传 距离 之 
间 不 存在 重 毒 区 域 (图 1) 。 采 用 NJ 法 构建 的 系统 
发 育 树 (图 2) ,每 个 物种 都 各 自 聚 为 一 文 , 文 持 率 在 
9996 以 上 ;结合 数据 库 中 已 公开 的 35 种 离 腹 寡 毛 实 
蝇 属实 晶 的 COI 基因 序列 ( 表 2) ,共计 47 种 实 蝇 ， 
构建 NJ 系统 发 育 树 , 每 种 实 晶 均 具 有 一 独立 分 支 
(图 3) 。 这 进而 证 实 ,DNA 条 形 码 技术 能 有 效 地 区 
A) MEE SOUS SCR SE TB 

近年 来 ,已 有 多 种 分 子 生物 学 技术 用 于 实 蝇 的 
检测 鉴定 。 如 扩 增 片段 长 度 多 态 性 技术 (amplified 
fragment length polymorphism, AFLP ) ( Kakouli- 
Duarte et al., 2001) RFLP 技术 ( 吴 佳 教 等 ，2004 ; 
Chua et al., 2010) , RAPD 技术 ( 张 亮 和 张 智 英 ， 
2007) 实时 荧光 定量 PCR 技术 ( 余 道 坚 , 2005) 、 基 
因 世 片 技术 ( 李 文 分 等 , 2008 、 微 卫星 标记 技术 
( microsatellite markers, SSR ) ( 姚 廷 山 等 ,2009 ) 、 核 
酸 测 序 技术 (Muraji and Nakahara, 2001; Smith et 
al., 2003; Zhang et al., 2010) 等 。 然 而 ,上述 检测 
鉴定 技术 各 有 优 缺 点 ,如 RFLP 技术 和 RAPD 技术 
具有 DNA 需求 量 少 多 态 性 丰富 、 对 基因 组 的 覆盖 
范围 比较 广 等 优点 ;但 是 , RFLP 技术 操作 繁琐 、 多 
态 性 检 出 率 低 ,而 RAPD 技术 的 稳定 性 较 差 、 对 实 
验 条 件 敏感 .可 信和 度 不 高 等 ,因此 限制 了 其 广泛 应 
Ho AE AFLP 技术 很 好 地 弥补 了 上 述 两 种 方法 的 
不 足 , 但 对 样本 DNA. 的 质量 要 求 较 高 , 且 所 需 设备 
的 费用 较为 昂贵 。 而 基因 忆 片 技术 虽然 可 以 满足 高 
量 的 需求 ,但 是 也 存在 技术 成 本 较 高 .操作 过 程 较 
为 复杂 等 问题 。 因 此 ,上 述 检测 方法 难以 满足 果蔬 
国际 贸易 及 国内 调运 过 程 中 快速 .准确 、 及 时 、 低 成 
本 的 检验 检疫 需求 。 尽 管 DNA 条 形 码 技术 已 用 于 
不 同 物种 的 分 类 人 研究 ,然而 ,该 技术 在 实 蝇 类 害虫 检 
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有 具 经 济 重 要 性 。 由 于 桔 小 实 蝇 复合 种 间 形 态 相 似 ， 
习性 相近 ,难以 快速 准确 地 进行 识别 鉴定 ,一 直 以 来 
都 是 分 类 学 家 研究 的 重点 。 本 研究 以 桔 小 实 蝇 、 杨 
桃 实 蝇 、 木 瓜 果 实 晶 .芒果 实 蝇 .菲律宾 果实 晶 .胡桃 
果实 晶 和 康 迪 果实 蝇 等 7 种 最 具 经 济 重要 性 的 桔 小 
实 蝇 复合 种 为 台 标 ,采用 DNA 条 形 码 技术 ,获得 3 
种 村 小 实则 复合 种 的 COI 基因 碱 基 序 列 ; 然 后 结合 
NCBI 中 已 公开 的 其 他 4 种 桔 小 实 蝇 复 合 种 的 COI 
基因 序列 ,以 邻接 法 构建 系统 发 育 树 (图 4) ,结果 显 
示 ,7 种 桔 小 实 晶 复 合 种 分 为 两 大 支 ,其 中 , 桔 小 实 
蝇 、 杨 桃 实 晶 、 木 瓜 果 实 晶 芒果 实 晶 和 菲律宾 果实 
蝇 聚 为 一 文 ,胡桃 果实 蝇 和 康 迪 果实 蝇 聚 为 另 一 文 。 
对 7 种 桔 小 实 量 复 合 种 的 所 有 COL 基因 序列 进行 种 
内 (5 种 ) 和 种 间 (7 种 ) 遗 传 距离 分 析 ,结果 显示 ,种 
内 遗传 距离 为 0.0006 ~ 0. 0080 ,平均 为 0.0059 , 其 
中 菲律宾 果实 晶 的 种 内 遗传 距离 最 小 为 0. 0006 , 杨 
桃 实 蝇 的 种 内 遗传 距离 最 大 为 0. 0078( 表 5 ) ;种 间 
遗传 距离 为 0. 0040 ~ 0. 0816 ,平均 为 0.0428, 其 中 
杨桃 实 蝇 与 胡桃 实 蝇 之 间 的 遗传 距离 最 大 为 
0. 0816 , 桔 小 实 蝇 与 菲律宾 果实 晶 之 间 .木瓜 实 蝇 与 
菲律宾 果实 量 之 间 的 遗传 距离 最 小 为 0. 0040 ( 
5)。 同 时 ,分 析 还 显示 桔 小 实 蝇 与 杨桃 实 晶 、 木 瓜 
实 蝇 、 芒 果实 蝇 、 非 律 宾 果 实 蝇 等 4 种 实 蝇 之 间 的 遗 
传 距 离 较 近 为 0. 0040 ~ 0. 0187 , 远 小 于 这 4 种 实 蝇 
分 别 与 胡桃 实 蝇 和 康 迪 果实 蝇 之 间 的 遗传 距离 
0. 0709 ~ 0. 0816 ; 而 胡桃 实 蝇 与 康 迪 果实 蝇 之 间 的 
遗传 距离 较 近 ,为 0.0165; 这 些 结果 说 明 , 桔 小 实 
Mà A MESE M ARK HR EAR SCRURT GEIE TE R SC 8 
具有 更 近 的 亲缘 关系 。 此 外 , 桔 小 实 晶 、 杨 桃 实 蝇 、 
木瓜 实 晶 菲律宾 果实 晶 等 4 种 实 晶 间 的 种 间 遗 传 
差异 小 于 条 形 码 有 效 性 的 检验 标准 0. 02 ( Herbert et 
al. , 2003) ,存在 种 间 差 异 比 种 内 差异 小 的 现象 ,本 
研究 结果 与 目前 国外 相关 研究 结果 一 致 (Mclnnis e 
al., 1999; Muraji and Nakahara, 2001), ix xk Hj] 
DNA 条 形 码 技术 不 能 有 效 地 区 分 桔 小 实 晶 复 合 种 。 
尽管 使 用 单一 一 种 分 子 标记 技术 对 该 复合 种 的 鉴别 
不 其 理想 ,但 通过 联合 应 用 DNA. 条 形 码 技术 、PCR- 
RLFP iX ,SYBR Green 实时 获 光 定量 PCR 技术 等 亦 





























































































































































































































测 识别 中 的 应 用 研究 尚 处 于 起 步 阶 段 , 如 Li 等 
(2011) 初 步 建立 了 包含 167 种 约 650 条 COI 基因 
序列 的 实 蝇 条 形 码 鉴定 系统 (Tephritid Barcode 
Identification System, TBIS) ,但 未 曾 进 行 系统 发 育 
分 析 ,而 且 该 系统 目前 还 不 宜 公 开 使 用 。 

桔 小 实 蝇 复 合 种 由 75 种 实 蝇 组 成 ,其 中 8 种 最 














可 有 效 地 区 分 部 分 桔 小 实 晶 复 合 种 ,如 桔 小 实 晶 、 杨 
桃 实 量 、 木 瓜 实 晶 .芒果 实 蝇 、 非 律 宾 果实 量 等 ( 余 
道 坚 , 2005; Chua et al., 2010) 。 然 而 ,目前 对 于 桔 
小 实 晶 复 合 种 中 村 小 实 晶 与 木瓜 实 蝇 杨桃 实 蝇 等 
种 类 的 界定 还 存在 争议 ,如 这 些 实 晶 在 形态 上 非常 
相似 ,而 且 存 在 生殖 隔离 不 完全 的 现象 ,不 同 种 类 的 
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雌雄 成 虫 可 互相 交配 并 产生 具有 繁殖 能 力 的 后 代 等 
( Melnnis et al., 1999) 。 

本 研究 采用 基于 COI 基因 序列 的 DNA. 条 形 码 
BUR VA E EXE SCRUR 21 种 具有 重要 经 济 意 
义 的 实 晶 为 靶 标 ,通过 种 内 、 种 间 遗 传 距离 计算 和 序 
列 聚 类 分 析 , 并 结合 NCBI 中 已 公开 的 部 分 实 晶 的 
COL 基因 序列 ,证实 了 DNA 条形码 技术 可 以 用 于 党 
见 实 蝇 类 害虫 的 识别 鉴定 ( 桔 小 实 蝇 复 合 种 除外 ) ， 
满足 口 赃 高 通 量 、 快 速 \ 准 确 、 及 时 鉴定 物种 的 需求 。 
尽管 本 研究 涉及 的 种 类 有 限 ,但 所 构建 的 实 蝇 类 
DNA 条 形 码 鉴定 技术 体系 ,为 有 效 阻截 潜在 实 蝇 类 
害虫 的 入 侵 和 局 部 发 生物 种 的 进一步 扩散 ,保护 农 
业 生 产 安 全 和 生态 安全 ,维护 国际 贸易 信誉 提供 了 
TUR SÉ 
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